Thermal Effects in the dynamics of disordered elastic systems 
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Abstract 

Many seemingly different macroscopic systems (magnets, ferroelectrics, CDW, vortices,..) can be described as generic 
disordered elastic systems. Understanding their static and dynamics thus poses challenging problems both from 
the point of view of fundamental physics and of practical applications. Despite important progress many questions 
remain open. In particular the temperature has drastic effects on the way these systems respond to an external force. 
We address here the important question of the thermal effect close to depinning, and whether these effects can be 
understood in the analogy with standard critical phenomena, analogy so useful to understand the zero temperature 
case. We show that close to the depinning force temperature leads to a rounding of the depinning transition and 
compute the corresponding exponent. In addition, using a novel algorithm it is possible to study precisely the behavior 
close to depinning, and to show that the commonly accepted analogy of the depinning with a critical phenomenon 
does not fully hold, since no divergent lengthscale exists in the steady state properties of the line below the depinning 
threshold. 
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1. Introduction 

Disordered elastic systems (DES) are ubiquitous 
in condensed matter. They are in particular rele- 
vant for magnetic [l|2|3l4j or ferroelectric [5|6j do- 
main walls, contact lines [7], fractures [819] , vortex 
lattices |10lllj , charge density waves [12] or Wigner 
crystals [13]. In these systems, one particularly im- 
portant question is to understand the response of 
the interface to an externally applied force, such re- 
sponse being directly measurable in all the above 
systems [14] . This problem is an extremely com- 
plicated one, since the competition between elas- 
ticity and pinning leads to the presence of many 
metastable states and thus to glassy properties. 

Despite this inherent complexity, we start now to 
understand well the dynamical properties of such 



systems at zero temperature. In particular, the sys- 
tem remains pinned up to a critical force Fc and 
starts sliding above this force. A particularly fruit- 
ful approach was to view this depinning transition 
in the light of critical phenomena [15], where the 
external force F would act as a control parameter 
analogous to the temperature leading to the transi- 
tion. The velocity v would be the analogous of the 
order parameter, finite on one side of the transi- 
tion and zero on the other side. This analogy imme- 
diately suggests the presence of critical exponents 
such as the ones defining the order parameter v oc 
(F — Fc)^ and the existence of a divergent length- 
scale ^ cx (F — Fc)~'^, that can be in that case iden- 
tified with the avalanche size when the system starts 
to move. These intuitive ideas are well confirmed in 
numerical and analytical studies of the T = sys- 



Preprint submitted to Elsevier 



26 September 2008 



terns. Given the success of the intuition gathered 
from the analogy with the critical phenomena, one 
is tempted to extend it to the case of finite tem- 
perature. In that case, thermal activation always al- 
lows to pass barriers ensuring a finite velocity at any 
force. Close to the depinning force Fc the finite tem- 
perature would thus be naively analogous to an ex- 
ternal field (such as a finite magnetic field) coupling 
to the order parameter (e.g. the magnetization in a 
magnetic system) and ensuring that this parameter 
is finite on both sides of the transition. This analogy 
suggests immediately two important consequences: 
a) there must exist a rounding of the depinning tran- 
sition due to the temperature and a thermal round- 
ing exponent V{Fc) ~ , as for standard critical 
rounding; b) a divergent lengthscale must exist on 
both sides of the transition, and at "short" length- 
scales the system must resemble the "critical" sys- 
tem, while above this divergent length the system 
must refiect the corresponding fixed point (high or 
low temperature for a normal critical phenomenon, 
zero force or high force in the case of the depinning 
transition). These consequences for the finite tem- 
perature behavior are quite striking and of course 
of direct experimental relevance, since experiments 
are always performed at a finite temperature. 

However, despite the success of the analogy with 
critical phenomena for T = 0, its extension to finite 
temperature is not guaranteed. It is thus important 
to check what remains of these properties. We ad- 
dress here these two points, and show that while a 
thermal rounding does indeed exist, there is no di- 
vergent lengthscale associated with the steady state 
properties of the line for F < Fc, showing that the 
analogy with a standard critical phenomenon does 
not hold and making it a challenge to develop an in- 
tuitive picture of the depinning transition at finite 
temperatures. 



2. Thermal rounding 

We first concentrate on the question of thermal 
rounding. 

Although analytical equations for the motion at 
finite temperature do exist [16j they are in practice 
very complicated to solve, and close to depinning no 
definite picture has emerged from them yet. However 
this question can be tackled by a numerical analysis 
of the equations of motion. We used a new technique 
allowing to avoid sample to sample fluctuations of 
the pinning force, which allowed us to obtain a very 
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Fig. 1. Steady state velocity against temperature for differ- 
ent system sizes. For large enough systems (L > 1024), all 
finite size effects are suppressed in the range of tempera- 
ture measured. This can be seen by comparing the data for 
L = 1024, 2048, that has been offset for clarity. The dashed 
line indicates the v ^ behavior, with the thermal expo- 
nent ~ 0.15 fitted using only L = 1024, 2048. After [IT]. 

good statistics and to extract the depinning expo- 
nent. Details and further reference to previous stud- 
ies of this problem in the literature can be found in 
Ref. [T7]. Fig. [T] shows the steady state velocity of a 
system at the critical force. 

Data shows a good scaling of the form v ~ T^, 
leading to the value ip ~ 0.15. More importantly 
one could show by measuring the structure factor of 
the interface that both the thermally-rounded and 
the r = depinning, display the same large-scale 
geometry, described by an identical divergence of a 
characteristic length ^ with the velocity of the form 
^ oc v~'^^^ , where i/ and /3 are respectively the T = 
correlation and depinning exponents. 

This remarkable fact fits very well universal scal- 
ing ideas from standard critical phenomena. It is 
important also to point out that so far there is no 
theoretical explanation of the value of the thermal 
rounding exponent (see Ref. TT] for a detailed dis- 
cussion on that point). Quite paradoxically the best 
estimate seems to be provided by a direct applica- 
tion of the formula for static critical phenomena, 
which of course has no justification in the case of 
such dynamical transition. From a more fundamen- 
tal point of view the scaling found above strongly 
suggests that the rounding of singularities that ap- 
pear in the renormalization study of the depinning 
[IB] are still controlled at finite temperatures by the 
velocity alone, as was the case for the T = depin- 
ning. 
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Fig. 2. The naive vision of the characteristic lengthscales and 
roughness of the line, in its moving frame. The analogy with 
critical phenomena would suggest a divergent lengthscale on 
both sides of the transition and a critical regime. This is in 
very different from what actually happens and is depicted 
in Fig. [3] 

The analysis of thermal rounding would thus com- 
fort the idea to extend the analogy with standard 
critical phenomena to finite temperatures. As we 
will see this is not the case when the characteristic 
lengthscales of the line are examined. 

3. Characteristic lengthscales 

A DES is characterized by a roughening exponent 
C that describes how the displacement from a flat 
configuration grow with the distance L between two 
points u{L) ^ . For a line in a two dimensional 
space, three roughness exponents can be easily de- 
fined. The equilibrium roughness exponent corre- 
sponding to zero applied force is C,eq = 2/3. If the 
line is moving very fast F ^ Fc '\t averages over dis- 
order and the system is thus equivalent to a pure 
thermal system. In one dimension, thermal fluctua- 
tions alone lead to a roughness exponent Qth — 1/2. 
Finally, at zero temperature, if the line is right at 
the critical force, it is characterized by the rough- 
ness exponent (^dep — 1-25 [18 . The analogy with 
a standard critical phenomenon would thus suggest 
that an observation of the roughness of the line (in 
the moving frame in which the line is at rest) would 
lead to a divergent lengthscale on both side of the 
transition separating the various roughness regimes 
as depicted on Fig. [21 However this picture disagreed 
with a functional renormalization group analysis of 
the small force (creep) regime, that was suggest- 
ing that (eg should be observed at short, not large 
lengthscales [IB] . 
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Fig. 3. Correlation length (squares) below the depinning 
transition, as computed by the algorithm described in the 
text. The line is a guide to the eye. Contrarily to what 
happens for F > Fc there is no divergent lengthscale showing 
up in the steady state properties of the line for F < Fc, in 
contrast with the naive picture that could be intuited from 
the analogy with a standard critical phenomenon. After |19) . 

Given the difficulty to analyze this question an- 
alytically or by using the standard numerical algo- 
rithms (Langevin dynamics for example) due to the 
extremely long computational times needed close 
to the depinning, we developed a new algorithm to 
tackle such a problem. Details on the algorithm can 
be found in fl9"20] . The basis idea is to use the fact 
that for very small temperatures, the system will 
spend most of its time in the configuration that has 
the largest escape barrier, since the escape time be- 
comes infinitely longer than for all the other config- 
urations. The steady state properties of the system 
are thus dominated by this particular configuration. 
This is the generalization for dynamics of the case 
of equilibrium thermodynamics where the partition 
function of the system is dominated at low temper- 
atures by the ground state of the system. Our algo- 
rithm allows by enumeration to find such a config- 
uration. Given the fact that the algorithm does not 
simulate the (extremely slow) dynamics close to the 
depinning, it can obtain quite quickly this configu- 
ration without suffering from the slowdown of the 
dynamics at low temperatures. 

The results are sketched in Fig. [3] They show 
clearly that no divergent lengthscale exists in the 
steady state properties of the line for F < Fc, in 
direct contrast with a naive extrapolation of our in- 
tuition on standard critical phenomena. In order to 
check that this result was no an artefact or our algo- 
rithm and that no subtle inversion of limits T 
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tion, we have also shown that the depinning transi- 
tion is not accompanied by a divergent lengthscale 
for F < Fc that would have bearings on the steady 
state properties of the line. This is in stark contrast 
with what happens in a normal critical phenomena, 
and directly prompts for a good description of the 
depinning transition. 
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Fig. 4. Characteristic length ^ as a function of the force / for 
different temperatures T increasing from the top to the bot- 
tom (T = 0.0025,0.005,0.01,0.025,0.05,0.075,0.1,0.15,0.2) 
as computed by Langevin dynamics. ^ grows monotonically 
with decreasing /. This is incompatible with the possibiUty 
of having a divergent steady-state length below threshold 
and is in full agreement with the results obtained by a dif- 
ferent algorithm and shown in Fig. |3] 

and F Fc could affect this result, we also spend 
considerable time recovering the crossover length by 
more conventional Langevin dynamics simulations. 
The results are shown in Fig.[?]and fully confirm the 
above results and the absence of divergent length- 
scale on the lower side of the transition. 

These results clearly shows that the intuition 
based on standard critical phenomena should be 
taken with a grain of salt and poses the problem of 
the theoretical description of the depinning transi- 
tion at finite temperatures. As for the case of the 
thermal rounding, although this information is in 
principle encoded in the FRG equations at finite 
temperature, extracting the physical behavior from 
these equations remain a very challenging task. 

4. Conclusion 

We have analyzed the effect of finite temperature 
on the depinning transition of a disordered elastic 
system. We have shown that at finite temperature 
the depinning transition is thermally rounded. The 
thermal rounding exponent ip :^ 0.15 needs a theo- 
retical explanation, since for the moment the only 
formula which seems to fit quantitatively is an un- 
justified extrapolation of the static critical rounding 
exponents for a standard critical phenomenon. Al- 
though the existence of thermal rounding would be 
in good agreement of the extrapolation of our intu- 
ition on critical phenomena to the depinning transi- 
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